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Constraints on /(R^kiR 1 ^ 1 ) gravity: An evidence against 
the covariant resolution of the Pioneer anomaly 
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We consider corrections in the form of AL(Rijti R ijkl ) to the Einstein-Hilbert La grangian. Then 
we compute the corrections to the Schwarszchild geometry due to the inclusion of this general term 
to the Lagrangian. We show that AL3 = ai{RijklR )^ gives rise to a constant anomalous 
acceleration for objects orbiting the Sun onward the Sun. This leads to the conclusion that ai = 

would have covariantly resolved the Pioneer anomaly if this value 



(13.91 ± 2.11) x 1CT 2B (— !— 

V ' \ meters 

of ai had not contradicted other observations 
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We notice that the experimental bounds on AL3 grows stronger in case we examine the defor- 
mation of the space-time geometry around objects lighter than the Sun. We therefore use the high 
precision measurements around the Earth (LAGEOS and LLR) and obtain a very strong constraint 



on the corrections in the form of AL(R ijk iR ljkl ) and in particular AL = a n (RijkiR ljkl ) n ■ This 



bound requires ai 
Pioneer anomaly. 



< 6.12 x 10 29 f — 7 — )s. Therefore it refutes the covariant resolution of the 

— meters ' 



I. DYNAMICS OF THE EMPTY SPACE-TIME 
GEOMETRY 

As we currently understand nature matters 1 deform 
their surrounding space-time geometry and gravity is a 
side-effect of this deformation. 

The simplest secular action capable of describing the 
deformation around matters' distribution, or equivalently 
the dynamics of the space-time, is the Einstein-Hilbert 
action: 



Seh = / d 4 x\/— det gR, 



(1) 



which happens to provide a very good phenomenological 
description of the gravitational phenomena within the 
solar system. Let us describe the solar system itself in 
a covariant language. In so doing we notice that the 
space-time geometry within the solar system possesses 
the following covariant characters: 

1. An almost vanishing Ricci scalar and tensor outside 
the world lines of the Sun and Planets, R = 0, 
R^v = 0. 

2. RavXvR^ > 1 



(meters) 4 



It should be highlighted that the first of the above as- 
signed properties to the Solar system is not exact due 
to the cosmological constant. The presence of a cos- 
mological constant of order A » 10~ 12 °t — - — seems 

'Planck 2 

to be the most economical description of the expan- 
sion of the universe which has been indicated by the 
distance-redshift relation. Ignoring the cosmological con- 
stant, however, seems to be a legitimate approximation 



when RaklR %3M > A 2 = 10" 



-104 



holds. Since 



Hjkl**' ^ meters 4 

RijkiR i: > kl > A 2 is violated beyond 10 7 AU from the Sun, 
a Ricci-flat geometry approximation to the geometry in 
the Solar system should be indeed very accurate. 

Note that the Newtonian gravitational interaction is 
the leading term in the effective gravitational interaction 
assigned to ([T|) . If we had been interested only in the ef- 
fective gravitational interaction then we could have sim- 
ply introduced some desired distance-dependent terms 
into the phenomenological effective gravitational poten- 
tial. "Distance" , however, is not a covariant quantity. 
We adhere to the standpoint that allows only generally 
covariant modification to ([I}. Once we add a modifica- 
tion, we can address what the modified term implies for 
the geometry around the Sun, Earth and even for time- 
dependent solutions. 2 This perhaps enables us to employ 
the high precision measurements around the Earth to test 
the validity of a covariant correction that is suggested by 
some accurate observations around the Sun. 

From a phenomenological standpoint, the Einstein- 
Hilbert action is only a model compatible with obser- 
vational quantities and experimental data. Unborn or 
approaching finer observations and preciser experiments 
may lead to some corrections to the Einstein-Hilbert ac- 
tion. These corrections can be any scalar constructed 
from the Riemann tensor and its covariant derivatives. 

From the theoretical standpoint, the simplest example 
of the corrections to the action we may consider is an 
arbitrary functional of the Ricci scalar and Ricci tensor: 



'Electronic address: exir@theory.ipm.ac.ir 
1 Non-vanishing energy momentum tensor. 



2 This standpoint also allows us to address what observed his- 
tory of the cosmological evolution of the universe implies on the 
time-independent solutions. This suggests that what causes the 
cosmological inflationary paradigm provides a lower bound on 
the mass of black holes [l[ . 
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L = R + eAL + 0{e 2 ), 
AL = f(R,Rij), 



(2a) 
(2b) 



where e encodes the perturbative expansion. The La- 
grangian given by ([2]) leads to the following equations of 
motion 



R^u + ekg llu f(R, Rij) + eO Ml ,[ 



df(R,IUj) df(R, Ri 



dR 



dR k 



0(e 2 ), (3) 



where is a linear operator acting on its arguments, 
and k is a constant number. We should solve ([5]) in a 
perturbative fashion: 



(4) 



where gj$ is a Ricci flat metric. Inserting (J4j) in ((3J leads 
to 

rt°Hg$) + kf(R,R ij )\ Ri ^ g$ + 



df(R,R l3 ) 
+°W JT5 lflii=o> 



df(R,Rjj) 
dR lrt,J - u ' dR M 



k,=o] = 0, (5) 



where is a linear second order differential operator. 
When f{R, Rij) has an expansion in term of its variable 
around Rij = then its partial derivative with respect 
to R or R^ either vanishes or diverges at Rij = 0. 3 If 
the partial derivatives vanish then a perturbative solution 
exists. For such case since is a linear operator then 
([5]) simplifies to 



□ ( %S) + fc/(o,o) 3 w = o. 



(6) 



Note that the above equation could have been obtained 
from the variation of the Einstein-Hilbert action in the 
presence of a tiny cosmological constant 



S = J d A x^J- det g(R + eA) + 0(e 2 ) 



(7) 



where A oc /(0, 0). The exact solutions of ([7]) are known. 4 
So the inclusion of the corrections in the form of a func- 
tional of the Ricci scalar and the Ricci tensor does not 
give rise to not-yet investigated perturbative corrections 
to the space-time geometry around the Sun. The same 
conclusion holds for any functional of the Ricci scalar, 
Ricci tensor and their covariant derivatives provided the 
functional has an expansion in terms of its variables 
around a Ricci flat geometry. 



3 Examples: The partial derivatives of f(R) = R 2 vanishes at 

R = while the partial derivative of f(R) = R^ diverges. 

4 Besides the Solar and cosmological constraints on f(R) or R n 
make them not attractive [2l| . Also look at the appendix A for 
a discussion on a misunderstanding of i gravity. 



A correction to the action thus would 'non-trivially' 
perturb the space-time geometry around a Ricci flat ge- 
ometry in case the correction involves the Riemann ten- 
sor per se. The simplest of such corrections is perhaps 
an arbitrary functional of the Riemann tensor squared: 



L = R + eAL + 0{e 2 ), 
AL = L(R m R^ kl ) , 

Examples of this form of correction include: 

AL l =axR i j H R ijk \ 
AL 2 = axiRijuR^ , 
AL 3 = ai(Ri jkl R ijkl )i . 



(8b) 

(9a) 
(9b) 
(9c) 



Having modified the dynamics of the space-time in small 
distances, the Heterotic and type I string theories leads to 
corrections to the Einstein-Hilbert action which includes 
ALi. 5 Note that ai is a dimensionless parameter, a 
constant number. Thus ai represents a possible struc- 
ture constant for the space-time geometry. On the other 
hand, the dimension of ai is negative with respect to the 
dimension of a±. AL3 thus may require a modification 
of the dynamics of the space-time in very small Riemann 
curvatures. 

This work aims to identify the best experimental lim- 
its and observational bounds on AL(RijkiR I -' kl ) and es- 
pecially AL — a n {RijkiR i: ' kl ) n regardless of what theory 
governs the dynamics of the space-time in a very high 
or low but not-yet achieved curvature. The work will be 
organized in the following order: 

In the second section, we will add a general correction 
in the form of AL(RijkiR^ kl ) to the Einstein-Hilbert La- 
grangian. After that we will compute the corrections to 
the Schwarszchild geometry due to the inclusion of this 
general term in the Lagrangian. 

In the third section, we will obtain the effective modi- 
fied Newtonian gravitational potential for space-craft in 
the spherical and static extrema of R + eAL + 0(e 2 ). We 
show that AL = ai{RijkiR K culminates in a con- 
stant anomalous acceleration for the objects orbiting the 
Sun toward the Sun. This leads us to the conclusion that 
qi = (13.91 ± 2.11) x 10 _26 ( ss ^ s )^ would have covari- 
antly resolved the Pioneer anomaly if this value of a 1 had 
not contradicted other observations. We notice that the 
experimental bounds for such correction grows stronger 
for the space-time geometry around objects lighter than 
the Sun. 

In the fourth section, we will utilize the high precision 
measurements around the earth (LAGEOS and LLR) to 
obtain a strong limit on the corrections in the form of 



AL 



Oi n {R-ijkiR l ^ kl ) n ■ This bound requires ai to be 



5 In the compactification of type II string theories, we also get 
non-perturbative world-sheet corrections in the form of L\. 
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smaller than 6.12 x 10 29 ( mc t crs ) ° , therefore, it clearly 
refutes the covariant resolution of the Pioneer anomaly. 

In the fifth section, we will note that AL2 = 
tti (-Riifei -H* ) 3 gives rise to an effective logarithmic 
gravitational potential. We will discuss whether AL 2 = 
ai{RijkiR l; ' kl ) 1 can describe the anomalous flat rota- 
tional velocity curves of the spiral galaxies. We show that 
a simple correction in the form of AL2 = ai (RijkiR 1 ^ 1 )? 
is not either in agreement with the high precision mea- 
surements around the Earth, or fails to describe the flat 
rotational velocity curves of the spiral galaxies. 

In the last section, we will provide a summary of the 
results. 



II. GENERIC BUT SIMPLEST CORRECTIONS 
AROUND THE SUN 



The Schwarszchild metric is an isotropic and static so- 
lution to the Einstein-Hilbert action. In four dimensions, 
in the standard preferred coordinates, it reads 



ds 2 — 



2 (1 



-)dt* + 
r 1 



— +r' i (dr + sin 0dr), (10) 



to obtain 



Eq l3 = r.h.s. of CD = Eq^ + eEqg + 0{e) , (13a) 



(i) 



Eq£> = R ij --g ij R, 



EqW=2R i a ^R^ d£ W' 



dx 



(13b) 



g l] C{TV) (13c) 



, dC[x\ 

dx '-= w2 



whereafter 1Z 2 = i? AI „ ??7 i? Miyl ' 7 is inferred. 

In the following lines from the outset, we consider the 
perturbations around the Sun in the standard preferred 
coordinate: 



ds 2 = -A(r) c dt 2 + B(r)dr 2 + r 2 {d6 2 + sin 2 6d<f) 2 ) . 



A(r) = (l- T f)(l + ea(r)) + 0(e 2 ) . 
B(r) = ^(l + eb(r)) + 0(e 2 ) . 

Now let us utilize 



TV 



+ ■ 



(14a) 
(14b) 

(14c) 



(15) 



and define 



where r/j = 2 ~2 in which G is the Newton constant, m 
represents the mass and c stands for the speed of light. 
We are interested in the space-time geometry around the 
Sun, we thus set m = M & = 1.98 x 10 30 kg for which 
r h = r — 3km. 

As argued in the first section, since the Einstein ten- 
sor vanishes for the Schwarszchild metric, the geometry 
around the Sun receives corrections in case the correc- 
tions to the action involve the Riemann tensor per se. 
We consider the simplest form of these corrections, the 
ones which are generic functional of the Riemann tensor's 
square: 



S = J d i x Llg^iR^r,-,] + 0(e 2 



(11a) 



£[<? M *,iW] = %/-det ff (iJ + eClR^BT™]), (lib) 

where e is the parameter of the expansion and £ is a 
generic functional. Computing the first. variation of pip 
with respect to the metric, we obtain [3] 



= 



8L 

dgij 



8L 



OR 



-R 



-2V Q V/3 



dL 



(12) 



where partial derivatives are taken assuming that and 
R^vrn are independent variables, and the partial deriva- 
tive coefficients appearing in (|12|) are uniquely fixed to 
have precisely the same tensor symmetries as the var- 
ied quantities. Note that our conventions are such that 
R<peipe as well as the Ricci curvature scalar are positive 
for the standard metric on the two- spher e. W e then in- 
sert the explicit form of L presented in (jllb[) into (TTJ 



L(r) = C[TZ 2 } 
L(r) ~ dC[x] 



dx 



724 



(16) 

L\r) + 0(e). (17) 



It follows that the non-vanishing and independent com- 
ponents of (fT3|) for (fT4|) are 



±E$> = ^(l - ^®)((r - r Q )(a' + b')-a + b)+ 0(e 2 ) . 



r(r - r ) 



(18a) 

{(r-Te)a' + a-b) + 0(e 2 ) , (18b) 



Eqf] = f ((r - r )a" + 2a' - b > + Q(e 2 ) 



(18c) 



and 



Eq 



(i) 



c 2 (r — r ) r 6 

+ rL' + 6L 



2r ( 2r(r-r Q )L"-(2r-3r Q )Z') 



0(e 2 ) 



12 r 

{r _ re)Eq m = _^'+6rL_2r s (-3r s + 2r) ., 



(19) 



12 



+ 0(e 2 ) 

Eq$ = -~(r(r-r )L"-(2r-3r )L') 



rL' + 6L 2 

r 

12 



+ 0(e 



(20) 



(21) 



where the dependence of the various functions on r is un- 
derstood. Then (|13[1 gives rise to a non-homogeneous second 
order differential equation for a(r) and a non-homogeneous 
first order differential equation for b(r). In accordance with 



4 



the preceding studies of string world-sheet corrections to var- 
ious black holes @, 0, Q] we demand that the corrections must 
not diverge on the possible event horizon at r = ih, provided 
the contribution of C(TZ 2 ) remains bounded on the horizon. 
We find out that this precondition is satisfied by the following 
solution 



b(r) 
a{r) 



L"(x) 



dx 

2(r- r ) J IQ i* 



-4r / dx^±, (22a) 

[-12 r| {xL'{x) + L(x)) + 

(22b) 

+8 r Q x 2 L(x) + x 6 L(x) + 2x 4 b(x)^ , 

The above boundary conditions also reproduce the cor- 
rect results around ordinary stars. In order to illustrate 
this claim, we note that we are solving the equations of 
motion in a perturbative fashion. We note that RijkiR 1 ^ 1 
inside the Sun is larger than RijkiR 1 ^ 1 outside the Sun. 
Since the perturbation that we are interested in grows 
as the Riemann curvature decreases (see the end of third 
section) then the perturbation holds valid inside the Sun. 
Let (01 represent the perturbative solution. Inserting this 
perturbative expansion in the equation of motion leads 
to a second order non-homogeneous differential equation 
for g$ . Let us seperately write the equations of motion 
of gjxv for inside and outside the matter's distribution: 



= F[g. 
□ (0) .^ t = F[g 



(0)i 
in \ ' 

(0)i 
out\ i 



(23) 
(24) 



where indices are understood but not written and F 
stands for the non- homogeneous part. Outside the star 
we have g^J t = ffSchwarazcMld- The consistency of the per- 
turbation requires that g^ remains bounded inside the 
star. Requiring the existence of the metric and its first 
derivatives on the surface of the star then provides the 
physical boundary condition for g^ t . This boundary con- 
dition, however, requires first solving the equations inside 
the star. Therefore a method capable of reproducing the 
physical boundary condition (which does not need solv- 
ing the equations inside the star) is appreciated . 

Since the equation for g^ t is linear, g^t f° r a star 
reads 



7 (1) 

y out 



C 2-5coi 



o (1) 

"Non-homogcncous 



(25) 



where g^ v and g^on are solutions of homogeneous part of 



that respectively diverges or converges when they are 
extrapolated toward the Schwarzschild radius associated 
to the central mass. c\ and C2 should be chosen such 
that the metric inside and outside the star match each 
other. In the Einstein-Hubert action, and for a star that 
is formed by collapse of dust, c\ and C2 can only depend 
on the total central mass: C\ = C\{M), c 2 = c 2 {M). 
On the other hand c\ and ci are numbers while M is 
parameter with a dimension. Therefore ci and c 2 must be 



independent of the mass of the star. Subsequently, once 
a physical criterion fixes c\ and C2 for a spherical mass 
distribution then they are fixed for every spherical time- 
independent distribution of mass. Since c% = 0, c% = 1 
for a large black hole then it holds c\ = 0, C2 = 1 for 
the Sun. We thus conclude that (|22p also describes the 
perturbation around the Sun. 

It is worth noting that though the precondition we em- 
ployed does not fix all the boundary conditions, it fixes 
the radial dependent part of a and b. The constant parts 
in a and b do not affect the force exerted onto a space- 
craft, the quantity we need in the next sections. In par- 
ticular for eC[K 2 ] = a n {R IJ , v xr l R' lvXn ) n ', the large radius 
behavior of a(r) and b(r) simplifies to 



f _r_\3 — 6n 

o(r) ~ -4(12)"n(l - n)a n r%- 4n r ° 



3 — 6n 



K r ) - -Gj _ 3n)a(r) . 



where n ^ \ is presumed. 6 For n = \ one finds 

air) ~ b(r) ~ -(12)3qi ln( — ) . 

2 r Q 



(26a) 
(26b) 

(27) 



For sake of simplicity, we will consider only the corrections in 
the form of e£[1l 2 ] = a n (R f , l/ x v R'" /X,1 ) n but the results can be 
directly generalized to a general case. We also note that the 
perturbations increase as we go further away from the Sun. 
We, however, shall address the regime of the validity of the 
perturbation in the fourth section after having independently 
obtained the experimental bounds on a„. 



III. EFFECTIVE POTENTIAL AND THE 
COVARIANT RESOLUTION OF THE PIONEER 
ANOMALY 

The effective potential for spacecrafts in the spherical and 
static geometry of (|14p is 



Vef f (r) 



I 2 



E 2 



+ 



2r 2 B{r) 2c 2 A(r)B{r) 2B(r) 



(28) 



where the equatorial plane is chosen to be orthogonal to 
the angular momentum, and E stands for the energy (per 
unit rest mass) and / represents the magnitude (per unit 
rest mass) of the angular momentum. 

Inserting the e perturbative expansion series of A(r) 
and B(r) into (gSJ) yields 



Veffir) = v}° ) f (r) + eV}j ) f (r) + 0(e 2 ) : 
GM e I 2 GM Q l 2 



(29a) 
(29b) 



6 It is not legitimate to extrapolate (I2(it toward far infinity since 
it has been obtained after expanding the equations for regions of 
the space-time which meets R^x-qR^^ > 7 (meters)" ' 
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where a constant term is tacitly recognized in (|29b[) . 
We notice that V^j-(r) is the effective potential for the 

Einstein-Hilbert action. Thus it is eV3l(r) which gives 
rise to the anomalous accelerations for spacecrafts de- 
ployed to explore the outer solar system. 

The Pioneers 10 and 11, spacecrafts deployed for ex- 
ploring the outer solar system are reported to have expe- 
rienced a constant anomalous acceleration of magnitude 
a p = (8.74 ± 1.33) x lCT 10 ^ in the direction toward the 
Sun at a distance of 20 — 70AU (Astronomical Units) 
from the Sun 7 [2J. What causes this anomaly might be 
on-board systematics, but the smoking gun has not been 
found yet |3j. Let us see if a correction in the form of 
AL = L(RijkiR^ kl ) may resolve the Pioneer anomaly. 

The effective potential culminating to the observed 
constant anomalous acceleration of the Pioneers is 



VEarth 



eV$(r) = 



(30) 



from now on we assume that (|30|) is valid for r E — 
70AU. 

The anomalous acceleration is obtained from analyz- 
ing the Doppler shift of the electromagnetic wave that the 
Pioneers had been sending assuming that the space-time 
geometry around the Sun coincides to the Schwarszchild 
geometry. Deviation from the Schwarschild geometry af- 
fects the gravitational red/blue shifts and subsequently 
alters the anomalous acceleration assigned to the Pio- 
neers [H, [HI El- 8 We first assume that the deviation 
from the Schwarszchild geometry needed to describe the 
Pioneer anomaly does not significantly affect the gravi- 
tational red/blue shifts of the signals sent by the space- 
crafts. We confirm this assumption after having identi- 
fied the deviation from the Schwarzschild geometry. 

If the Pioneer effective potential (|30[) is due to gravi- 
tational effects of the Sun then (|29cp and ([ST)]) results 



E — c , , . c 2 , , 



(31) 



One notices that the Pioneers have classical velocity and 

E 2 ^i C - j « c 2 whereafter v stands for the radial veloc- 
ity of the spacecrafts with respect to the Sun. Therefore 
pip can be further approximated to 



2a 



ea(r) = -^r + O(^) . 



(32) 



where rearranging the terms is induced. Note that in 
these distances it holds ea(r) <C and the approxima- 
tion in (|32[) is much lesser than the error bar in a p . 

Eq. (I32[) identifies the deviation from the 
Schwarzschild geometry. In the following we are going 
to argue that this deviation does not significantly alters 



LOT 



Vpioneer 

> 



V4 



% - 


r 









FIG. 1: The acceleration of the Pioneer spacecrafts are iden- 
tified by analyzing the Doppler shift of the spacecrafts ' emis- 
sion. The gravitational red/blue shifts are part of this ana- 
lyze. Altering the geometry affects the gravitational red/blue 
shifts. Therefore, the back reaction of the change of geometry 
to the constant anomalous acceleration should be addressed if 
the anomaly is due to a covariant correction. 



the gravitational Doppler shifts to which the anomalous 
acceleration is assigned. The Pioneer spacecraft 9 sends 
a photon at time t in x = rp; onecr of frequency v\ in the 
local frame that is fixed at rpj oneer . The detector on the 
Earth recieves the photon on time t+T(t) with frequency 
of V2 in the local from that is fixed at x — rbetector > see 
fig. IIIII T{t) can be identified in terms of the initial 
position of the spacecraft, Earth and the detector. 

In the absence of gravity we assign energy of hv\ to 
a photon related to light's wave of frequency v\. In the 
presence of the week gravity, therefore, we should assign 
energy of hv + V e ff-^ to the total energy for the photon 
of frequency Hv where V e ff is the effective gravitational 
potential at where photon is localized. Requiring the 
conservation of the energy for the photon sent by the 
Pioneer yields 

v\ (1 + ^V e //[r Dctoctor(t+T) ]) = i/ 2 (l + ^■V e //[rp ioncer{t) ]) 

— ^(14//[rpiojieer(t)] — Veff [?"Detector(t+T)]) -(33) 



v-2 



Employing (|29[) then yields 
2 Av 



K// [detector (t + T)] - V§) [r Pioneer (t)] + 
+ ^{a[r Dete _ctor(t + T)] - a[r Pi „ ee r {t)]) . (34) 

We note that the time that a photon needs to travel from 
the spacecraft to the detector also has an e expansion: 

T = T (0) +eT (1) +0(e 2 ), (35) 

r D .{t + T) = r D {t + T il)) ) + er D {t + T {0) )T W +C>(e 2 ), 

where rr>. = ^Detector- Inserting (|3"5|) in [33] leads to 



Av ,Av 



dr D .{t) „ w 
C 2 dr |r=rD dt 



i' 

£0„ 



= (-)K-H,+.^| f=n ,.^T W + 



+ -f [rpioneer (t) - r Dctcctor (t + T (0) )] , (36) 

c z 



7 The major part in the error bar of a p is systematic. The statis- 
tical error is only 0.01 X lO -10 -^. 

8 I appreciates the comment made by the referee of the CQG that 
leads to adding a discussion on this issue. 



The Deep Space Network antennas tracked the pioneer space- 
crafts with a S-band signal at about 2.11 GHZ. The tracking 
was done by sending a signal from the Earth which the Pioneers 
were replying. The simple one-way Doppler shift analyze of this 
section, however, suffices for our conclusions. 
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where (4t)e.-h. is the prediction of the Einstein-Hilbert 
action. This states that ignoring the deviation from the 
Einstein-Hilbert action leads to a systematic error in ^ 
given by 



c 2 9r 



at 



T (1) + (37) 



+ ^[rpionc C r(i) — r D ctoctor(i + T (0) )] 



Since we are interested in a motion of a spacecraft, we 
consider only (|40|) . Then comparing (j4Tjl to (|32f identi- 
fies a to 



(13.91 ± 2.11) x 10~ 26 (- 



:)• 



(42) 

3 meter 

which would have covariantly resolved the Pioneer 
anomaly if it had not been in contradiction with the other 
observations. 



which in turn results to a systematic error in determining 
the acceleration of the spacecraft 



missing — C "77 ( )mts ~ . r» <i 



(38) 



We notice that w 



should be propor- 



tional to a„. Therefore 



This helps us 



. _ -lp . X C1U1C 1 ~ Lip A _ _ 

to obtain the order of magnitude of the terms present in 



(a) r . 



-REarthtt'Earth ^ 

(39) 

(8.74 ± 1.33) x 



Therefore, due to the error bar in a p = 
10 _10 f , it is legitimate to neglect a missm g. 

We note that the covariant resolution of the Pioneer 
anomaly gives rise to periodic terms (with periodicity of 
one day and one year) in the Doppler shift (|36|) . [l2j 
and [13| report that the residual of the fit with constant 
anomalous acceleration contains clear periodic terms. 
Ref. [IH argues that these periodic terms should be as- 
signed to the Earth and its atmosphere while [l2j dis- 
cusses that they are somehow fingerprints of what causes 
the Pioneer anomaly. The periodic terms reported in 

OBI are at order 01 4r ~ lift = 10 ~ 12 - The term 
of periodicity of one year in (|3"6")) is at order of magnitude 

^r\ yr « ^-2AU = 1CT 15 , and the term of periodicity 
of one day in (J36]) at order ^-\ Day « |2 x 6400fcm = 
10 -19 . Therefore no covariant resolution of the constant 
pioneer anomaly is able to account for the residual pe- 
riodicity of the fit with constant anomalous acceleration 
(note that this statement was derived only from Q32p). 

Now let us come back to the main issue of this section: 
how does (|3"2"j) help us to identify the correction to the 
action? We see that (|32|) beside (|22j) leads to an inte- 
gral equation for the correction to the Einstein-Hilbert 
action. Variation of this integral equation with respect 
to r leads to a non-homogeneous linear third order dif- 
ferential equation for L(r). Eq. (I26| then shows that one 
solution of this differential equation is 



Ol 
3 



X (i? M ^r J7 i?' J ^' 77 ) 3 



which leads to 



eai i 

3 



.12,1 8r 
-ai ( — )3 — 
3% ; 9 



6j(r) 



3 / ^ 
2 3 V ; 



(40) 



(41) 



It is worth noting that (|4"0|) is only one solution to the 
corresponding non-homogeneous linear third order differ- 
ential equation for L(r). Other solutions differ with (j40|) 
by terms which do not affect the motion of a spacecraft. 



IV. CONSTRAINTS FROM EARTH AND 
MOON 

In the previous section we have examined a general 
family of the covariant corrections and found the covari- 
ant correction capable of describing the Pioneer anomaly. 
We note that the corrections to the space-time geome- 
try given in eq. (|41j) would increase if the mass of the 



Sun decreases (recall that r Q 



2GM F 



). It implies that 



the accurate measurements of the geometry around the 
Earth would provide a strong constraint on the covariant 
corrections in the form of: 



S = j d 4 x sj- det g(R + Q, l (i? M ^ 7 i?"" 



(43) 



Recalling that (|43[) and (|29c|l hold a perturbative expan- 
sion like to that of (HU) and (|26H27|) around the Earth, 
we find that a probe with a classical velocity experiences 
an anomalous acceleration of magnitude a e (r), given by: 



I \ 2-6 

a e (r) = —a n r 
a n = 2(12) n n(l - n)c 



(44a) 
a n , (44b) 



at the distance r from the center of the Earth toward 
the center of the Earth in case that eq. (|43|) governs the 
dynamics of the space-time. To put it another, a satellite 
in a circular orbit experiences the following gravitational 
field around the earth 

GAfEarth 



Fa = 



+ a n r 



(45) 



where Msarth is the inertial (effective) mass of the Earth. 
The effective gravitational mass of the Earth is defined 
by 



GM, 



eff 
Earth 



(r) 



r 2 F< 



G 



GM Ea rth + a n r 4 ~ 



6n 



(46) 



In contrast to the Newtonian dynamic, the effective grav- 
itational mass is not a radius-independent quantity. An- 
alyzing the circular orbit of any satellite or the Moon 
identifies the effective gravitational mass of the Earth 
within that orbit. 

The accurate value of the mass ratio of the 
Sun/(Earth+Moon) from the Lunar Laser Ranging can 
be combined with the Solar GM and the lunar GM from 
lunar orbiting spacecrafts [l6| to give the effective gravi- 
tational mass of the Earth in an Earth-centered reference 
frame with the precision of one part in 10 8 : 



G-ME art h (dEarth-Moon) 



398600.443 ± 0.004- 



, km 3 



(47) 



7 




/ "Break \ 

g 



0.1 0.2 0.3 0.4 0.5 0.6 




FIG. 2: The blue region is the range of parameters describing 
the covariant corrections to gravity which are consistent with 
the accurate measurements of the space-time geometry around 
the Earth (|50|) . The dots represents the values which will 
be needed to covariantly resolves the Pioneer anomaly or the 
anomalous velocity curves of the spiral galaxies. The thickness 
of the dots corresponds to their uncertainties. 



where ^Earth-Moon is the distance between the Moon and 
Earth [l7[ . The effective gravitational mass of the Earth 
has also been measured by various artificial Earth satel- 
lites [lH, including the accurate tracking of the LAGEOS 
satellites orbiting the Earth in nearly circular orbits with 
semimajor axes about twice the radius of the Earth: 



ri, , LAGEOS \ 
LfJWEarth (2r E arthJ 



= 398600.4419 ±0.0002 - 



kiri 



(48) 



where rEarth stands for the radius of the Earth [19(. We 
note that within the errors (f48|) is compatible with (|47|) . 
Recalling (|46|) beside comparing (|47|) to ([48]) . therefore, 
leads to 

|a»| ((d Earth -Moon) 4 - 6n - (2r Earth ) 4 - 6 ") < 0.004^ . (49) 



where \a n \ stands for the absolute value of a r , 
(|49| can be approximated to 



For n < | , 



|Q.n|(dEarth-Moon) 4 6 " < 0.004 

S 



-r (50) 



Using (|44bp . ([50]) results 



< 



0.004 



km 3 I GM Ea 



\l-2n 



2(12)"n(l - n)c 2 (d E arth-Moo„) 4 - 6n 



(51) 



The values of |a„| which meet ([ST]) are illustrated in 
FigH for 0.05 < n < 0.65. The limit on |ai| is 

lail < 6.12 x 10- 29 ( — i — )§ . Therefore a p t = (13.91 ± 

1 t 1 — ^ meters > 4 ^ 



FIG. 3: We have used the bounds on a n which is presented in 
Fig. 1 to find how far from the Sun the perturbation remains 
valid. The perturbation in terms of a„(i? I - ,fc ' Rijkl) is valid 
below the blue continuous line. Below the dashed line, the 
perturbation around the Schwarszchild geometry is valid. Note 
that the y-axis is both logarithmic and represents \og lg ( r A r ff k ) 
where r Break is the minimum distance from the Sun that the 
'corresponding' perturbation breaks in. 



is in agreement with other independent studies: the pre- 
cession of the longitudes of perihelia of the solar planets 
[lij or the trajectories of long period comets [l5[ have 
not been reported to experience an anomalous gravita- 
tional field toward the Sun of the magnitude capable of 
describing the Pioneer anomaly. 

Having obtained the experimental bounds on a n (|51D . 
we would like to find the minimum distance from the 
Sun that the perturbation breaks in. 10 Note that we have 
assumed that the space-time geometry has a perturbation 
around the Schwarzschild geometry, 

- g tt = A(r) = (1 - !® )(1 + ea(r) + 0(e 2 )) , (52) 
r 

wherein we have assumed that the Schwarszchild geom- 
etry describes the space-time geometry with a very good 
approximation. The existence of the perturbation means 
that ea(r) <C 1. The Schwarzschild geometry remains 
a good approximation if ea(r) <C The perturbation 
breaks when ea(r) w — or ea(r) w 1. Using the bounds 
on a n we then obtain 



ea(r) 



r 1 + 1 12»4n(n - 1) a n r%- 4n U 



3 — 6n 



1 



ea(r) ~ > — « ( ; ; s — r— 

r r Q A 12"4n(n - 1) a n r 2 Q ~ 4n 



(53) 



2.11) x 10" 26 ( 



which is needed to covariantly re- wne re a„ is such that the bound in (51]) is saturated, and 



solve the Pioneer anomaly, is clearly not compatible with 
the accurate measurements around the Earth. So the Pi- 
oneer anomaly can not be covariantly resolved within the 
general family of corrections we have considered. This 
supports the idea that the Pioneer anomaly is on board 
systematic or due to non-gravitational effects. This idea 



in the last line it is assumed that n ^ Fig 2. plots 



10 I thank the comment of the referee of CQG that leads to adding 
the subsequent paragraphs in this section. 
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f. The first observation is that the Perturbation remains 
valid inside the Solar system. It is also interesting that 
the perturbation around the Schwarszchild metric breaks 
before perturbation in a n (RijkiR l -' kl ) n . This means that 
there exist some regions where in R + ea n {RijklK ) n + 
0(e 2 ) is perturbative in the sense that terms of order e 2 
can be consistently neglected while the dominant term in 
R + ta n {R ljkl R l3kl ) n is a n (R l3k iR l ° kl ) n not R. It might 
be interesting to solve the "exact" equations in these re- 
gions. 

We would like to reemphasize that the combined LLR 
and LAGEOS measurements provide indeed a strong con- 
straint on the form of the covariant correction. In order 
to further illustrate the power of this constraint let us 
investigate if it is satisfied by the covariant corrections 
proposed in ref. [33J. Ref. [331] considers a family of f(R) 
gravity in the presence of the cosmological constant be- 
fore showing that in the vicinity of the Sun, there exists 
a set of f(R) corrections capable of the describing the Pi- 
oneer anomaly. It shows that the effective gravitational 
acceleration is 

GMq 

^gravity ^2 ^constant 7 v^* / 

where the second term is "a constant acceleration [while 
is] independent of the [considered central] mass" . It then 
sets a constant = a p = (8.73 ± 1.33) x lCr 10 f . Ref. [H 
however has not considered the implication of the covari- 
ant correction they studied to the space-time geometry 
around the Earth. ([54]) implies the following effective 
gravitational acceleration around the Earth: 



^gravity 



GM Ea 



■Hi 



(55) 



which leads to the following effective gravitational mass 
of the Earth 



Gi^l = GM Ea 



uarth + a p r . (56) 

The combined LAGEOS and LLR measurements then 
requires 

, km 3 



ap(dEarth-Moon) 2 < 0.004- 



(57) 

while a p (d E arth-Moo„) 2 = (0.129 ± 0.020)^. Therefore 
the Earth-Moon system also refutes the covariant resolu- 
tion of Ref. [33[ for the Pioneer anomaly. 

A lesson we should learn here is that "phenomenolog- 
ically good" covariant corrections to the action which 
remain perturbative at the 'close' vicinity of a spheri- 
cal central mass seem to be those whose predicted cor- 
rections to the space-time geometry had decreased if we 
would have decreased the central mass. The perturbative 
studies of f(R i jkiR l -' kl ) [this paper] and f(R) [33j show 
that this criterion is not satisfied in general. It is inter- 
esting to systematically study what kind of corrections 
meets this criterion. In so doing, perhaps, it is inter- 

esting to study /( v R %3 * 1 v i R ]ki m ) that is suggested by 
the covariant resolution of the anomalous flat rotational 



V. ON THE ANOMALOUS ROTATIONAL 
VELOCITY CURVES OF THE SPIRAL 
GALAXIES 



Eq. (1271) an d (|28|) demonstrate that a correction in 
the form of AL = a^RijkiR 1 ^ 1 )? leads to the following 
effective perturbative gravitational potential 



Veff(r) = 



GM a 



- 2y/?>ai c 2 ln( — ) + 0{a\), (58) 



around any spherical static distribution of matter of 



total inertial mass M g 
ikt 



suggests that AL 



&i(RijkiR t:> ) 2 might have a chance to resolve the 
anomalous flat rotational velocity curves of the spiral 
galaxies without considering dark matter. Mathemat- 
ically speaking, the current work does not approve or 
reject this suggestion, due to the following reasons: 

1. The anomalous velocity curves of the spiral galaxies 
occur at the boundaries of the spiral galaxies. The 
matter's distribution inside the presumed galaxy 
is disk-like and not spherical for the stars in the 
boundary of the spiral galaxies. ([2"T]) and ([2"5]) are 
derived for a spherical distribution of matters. 

2. The anomalous velocity curves of the spiral galax- 
ies are not small deviation from what Newtonian 
gravity predicts. The exact solutions of the mod- 
ified action might precede a possible resolution of 
the anomalous velocity curve. 

Despite the above obstacles we tend to examine whether 
a value of ai compatible with (|50[) has a chance to 
describe the flat rotational velocity curves of the spi- 
ral galaxies. In so doing, let us extrapolate ([58]) to- 
ward the boundary of a typical spiral galaxy of mass 
AI g = 10 12 M Q . This generalization leads to the following 
relation for the velocity of the stars moving on a circular 
orbit around the center of the galaxy 



GAL 



+ 2%/3aic 2 + 0(a!), 



(59) 



where v stands for the velocity of the star around the 
center of the presumed galaxy. Examining the rota- 
tional curves of the spiral galaxies FigfJ] we see that 
the constant asymptotic velocity can be approximated 
by 200 ± 50^ in large scales. In these distances, the 
first term of ([59]) is small, thus ([59]) implies 



(60) 



a\ = (13.61 ±6.39) x 10" 



is needed to describe the constant velocity of the within 
the borders of a typical spiral galaxy. The high precision 
measurement around the Earth (]50[) . however, requires 



(61) 



a 3 , < 6.67 x 10" 



curves of the spiral galaxies [3J 



therefore, discards (|60|) and implies that AL — 
(RijkiR 1 ^ 1 )? has no chance to describe the flat rotational 
velocity curves of the spiral galaxies. 
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FIG. 4: Rotational curves of spiral galaxies obtained by com- 
bining CO data for the central regions, optical for disks, and 
HI for outer disk and halo J2(£j. 



It is worth noting that the Riemann scalar curvature 
in the Solar system and around the Earth satisfies 



7.53 x 10' 
meter 4 



^ Riwri^R 



(62) 



while in the regime where the anomalous rotational cur- 
vature of the spiral galaxies happens, it satisfies 



10- 



meter 4 



10- 



meter 4 



(63) 



where a galaxy with a central mass at order 10 12 Af Q 
with the boundary of about lOOkpc is alleged. If we as- 
sume that a simple action in the form of a polynomial 
in terms of the Riemann tensor dictates the dynamics of 
the space-time in both of the regimes given by (|62[) and 
([63|) we observe that the value of ai will be needed to 
describe the anomalous flat rotational curve of the spiral 
galaxies (|60[) is not in agreement with the accurate mea- 
surements of the space-time geometry around the Earth, 
as illustrated in Fig|2] as well. We, however, lack ex- 
perimental justification or observational data supporting 
this assumption. Any functional of the Riemann tensor 
squared , Q[K 2 ] = Q[R^ VV1 R^], which becomes suf- 
ficiently small for (f6"2")) but constant for can be uti- 
lized to suggest the following phenomenological action for 
gravity: 



S = 



j d 4 x \f— det g(R + a' 



Q[H 2 ] (iW-R^" 7 )^), (64) 



which has a chance not only to be consistent with the So- 
lar system's data but meets (|60p as well. The exact solu- 
tions of (|60|) precedes reaching a concrete conclusion on 
the validity of the above suggestion. Addressing the exact 
solution of (|60|) or a similar action in which i? AI „ ??7 i? A " yl ' 7 
is replaced with the Gauss-Bonnet Lagrangian lays out- 
side the scope of the current work. 



to the Einstein action would 'non-trivially' perturbs the 
space-time geometry around a Ricci flat geometry in case 
the correction involves the Riemann tensor per se. After 
that we have considered the simplest family of these cor- 
rections; the corrections which are arbitrary functional 
of the Riemann tensor's squared: 



L = R + eAL + 0(e 2 ), 
AL = L(R ijkl RV kl ), 



(65a) 
(65b) 



Then we have computed the corrections to the 
Schwarszchild black-hole in an asymptotically flat 4D ge- 
ometry for a general AL. 

We have observed that AL3 = a^{RijkiR gives 
rise to a constant anomalous acceleration for objects or- 
biting the Sun onward the Sun. This leaded us to the 
conclusion that ai = (13.91 ± 2.11) x lO -2 ^^^)^ 
would have covariantly resolved the Pioneer anomaly if 

this value of ai had not contradicted with other obser- 

3 

vat ions. 

We have shown that the experimental bounds on AL3 
becomes stronger in case we examine the deformation of 
the space-time geometry around objects lighter than the 
Sun. We therefore have used the high precision mea- 
surements around the Earth (LAGEOS and Lunar Laser 
Ranging) and have obtained a strong constraint on the 
corrections in the form of AL(R i jf t iR l: ' kl ) and in par- 
ticular AL — ctniRijkiR 1 ^ 1 )™ ■ R is interesting that 
the high precision measurements around the Earth pro- 
vide a strong constraint on the possible correction to the 
Einstein-Hubert action. 

The high precision measurements around the Earth re- 
quires ai < 6.12 x 10~ 29 ( mc t crs ) 5 ; therefore, they refute 
the covariant resolution of the Pioneer anomaly. So the 
Pioneer anomaly can not be covariantly resolved within 
the general family of corrections we have considered. 
This supports the idea that the Pioneer anomaly is on 
board systematic or due to non-gravitational effects. 

We also have noted that AL 2 = ai(RijkiR^ kl )" s gives 
rise to an effective logarithmic gravitational potential. 
We have raised the question if AL 2 = ai(RijkiR i: ' kl )^ 
may be useful in describing the anomalous flat rota- 
tional velocity curves of the spiral galaxies, before having 
proved that a simple correction in the form of AL 2 — 
a^iRijkiR 1 ^ 1 ) 1 is not either in agreement with the high 
precision measurements around the Earth, or can not 
describe the flat rotational velocity curves of the spiral 
galaxies. 
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APPENDIX A: SOLAR SYSTEM TESTS DO NOT 
RULE OUT 1/R GRAVITY 

There exists a debate in the literature on the consis- 
tency of 11 



S = 



1 



16ttG 



for n 

example look at 
In particular ref. 



J d 4 xy/-detg{R-^) + J d 4 xy/-det gL m 

(Al) 

10 26 meters with the Solar 

H IS H HI EE [27 




claims that properly matching the 



metric inside and outside the Sun rules out gravity. 
A sharp inspection of [22j], however, reveals that [22j has 
not properly studied the equations. This was also noticed 
by [26] . In the following we would like to clarify why the 
conclusion of (22j is not right in addition to demonstrat- 
ing the source of this wrong conclusion. 

Contracting the equations of motion of (|Alj) with the 
inverse of the metric yields 



10~ 31 meters 2 . Therefore even the non-homogeneouty of 
the matter's distribution in the Sun does not produce a 
significant deviation from x = 1 inside the Sun. Recalling 
that the Einstein-Hilbert gravity holds x — 1 besides 
extraordinarily small deviation from x = 1, we conclude 
that the Einstein-Hilbert action quite-perfectly describes 
the physics in and outside the Sun for (jAljl . 

Now let us inspect what leads [22( to the wrong con- 
clusion. Ref. [22] defines a function by 



R 2 



(A6) 



which we refer to as the C-function. The authors then 
take granted that the C-function encodes the deviation 
from the vacuum solution even outside the matter's dis- 
tribution. This means that the authors fail to realize 
that the C-function can be identically zero outside the 
matter's distribution. In other words the deviation from 
the vacuum solution might be encoded in other scalars 
rather the Ricci scalar or equivalently the C-function. 
They than rewrite (|A2|) in terms of the C-function 



□c- 



jjL 2 C 



8ttGT 



(A7) 



R 2 3 



^ _ 8irGT 
~R ~ 3 



(A2) 



where T — g^T^ and the speed of light is set one. Now 
let us define a new variable, x, through R = —8irGx. 
Rewriting (|A2|) in terms of x and rearranging the terms 
yields 



x = 1 



3a* 4 



inGT)' 



K4 + 



1 



8ttGT 



□-^+0(VT)). (A3) 



Note that the Einstein-Hilbert gravity holds x = 1. Wo 
can obtain the order of magnitude of T for the Sun by 
T « 4,fff1 where Rq is the radius of the Sun. Knowing 

the order of magnitude of T, we can obtain the order of 
magnitude for deviation from x = 1 in ()A3|) : 



3/ 



3H 4 



^GTf 



\ R| (speed of light) 2 > 



3r| 



(A4) 



where 



2GM 



is the Schwarszchild radius 



(speed of light) 2 

associated to the mass of the Sun. Using Rq — 1.39 x 
10 9 meters, r Q = 3km and /j w 10~ 2e meters, the order 
of magnitude of the deviation from x = 1 in (|A3j) reads 

Note that the order of magnitude of the coffiecient in the 
front of □ A- and similar terms in 11 lo — v — 



1S 7^ X (8 ;:(■■/ 'i- 



before approximating it to 



V 2 c- 



8ttGT 



(A8) 



It then appears that the authors assume that the C- 
function and its derivatives are continuous on the surface 
of the Sun. But in the Enistcin-Hilbcrt gravity what re- 
mains continuous on the boundaries are the metric and 
its first derivatives. For example we know that the Ricci 
scalar or equivalently the C-function is not continuous on 
the boundary. Therefore instead of choosing C-function 
as what [22j has chosen, we must choose it in the follow- 
ing way 

• Outside the Sun (r > Rq), c = 0. Note that c = 
solves the equation outside the Sun. 

• Inside the Sun we must find a solution of (|A7j) that 
remains bounded inside the star. 

The above choice means that within the electrostatic ap- 
proximation to the equations - (|A8[) -. the surface of the 
Sun effectively plays the role of a conducting surface ac- 
commodating some amount of 'charge' that completely 
cloaks the 'charge' inside the Sun. Setting c = outside 
the Sun leads to the Schawrszchild metric in the Solar 
system which is in agreement with observation. There- 
fore 1/R gravity is not ruled out at least due to reason 
addressed in 12211. 
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11 This appendix is added as part of the response to a criticism 
made by the referee of CQG. 



[1] Q. Exirifard, Lower bound on the mass of a black hole, 
larXiv:0807.' 3246 [gr-qc]. 

[2] M. M. Nieto and S. G. Turyshev, Finding the origin of [18] 
the pioneer anomaly, Class. Quant. Grav. 21 (2004) 4005, 
|gr-qc/0308017] . [19] 

[3] J. D. Anderson et al.,The search for a standard explana- 
tion of the pioneer anomaly, Mod. Phys. Lett. A17 (2002) 
875, |gr-qc/0107022] . [20] 

[4] V. Iyer and R. M. Wald, Some properties of Noether 

charge and a proposal for dynamical black hole entropy, [21] 
Phys. Rev. D50 (1994) 846 , [gr-qc/9403028] . 

[5] C. G. Callan et al, , R. C. Myers, and M. Perry, Black [22] 
holes in string theory, Nucl. Phys. B311 (1988) 673. 

[6] Q. Exirifard and M. O'Loughlin, Two and three 

loop alpha-prime corrections to T-duality: Kasner and [23] 
Schwarzschild, JHEP 12 (2004) 023 , [hep-th/0408200] . 

[7] Q. Exirifard, The world-sheet corrections to dyons in the [24! 
Heterotic theory, JHEP 03 (2007) 074, [hep-th/0607094| . 

[8] M. Blau, Lecture notes on general relativity. 

[9] M. T. Jaekel et al. , Gravity tests in the system and the 

Pioneer anomaly, Mod. Phys. Lett. A20 (2005) 1047; [25] 
Post-Emsteinian tests of linearized gravitation, Class. 
Quantum Grav. 22 (2005) 2135. 
[10] M. T. Jaekel and S. Reynaud, Post-Einsteinian tests [26 
of gravitation, Class. Quant. Grav. 23 (2006) 777 
[arXivigr qc/0510068j . [27] 
[11] M. T. Jaekel and S. Reynaud, Radar ranging and 

Doppler tracking in post-Einsteinian metric theories [28 
of gravity, Class. Q uant. Grav. 23, 7561 (2006) 
[arXiv:gr-qc/0610155] . 
[12] A. Levy, B. Christophe, P. Berio, G. Metris, J. M. Courty 

and S. Reynaud, Pioneer Doppler data analysis: study of [29 
periodic anomalies, arXiv:0809.2682 [gr-qc]. 
[13] J. D. Anderson, P. A. Laing, E. L. Lau, A. S. Liu, 

M. M. Nieto and S. G. Turyshev, Study of the anoma- [30] 
lous acceleration of Pioneer 10 and 11, Phys. Rev. D 65 
(2002) 082004 [arXiv:gr-qc/0104064| . 
[14] Lorenzo Iorio, Can the Pioneer anomaly be of gravita- [31 
tional origin? A phenomenological answer, Found. Phys. 
37 (2007) 897. 

[15] M. M. Nieto, S. G. Turyshev and J. D. Anderson, Di- 
rectly Measured Limit on the Interplanetary Matter Den- [32 
sity from Pioneer 10 and 11, Phys. Lett. B 613 (2005) 
11. [33" 

[16] A. J. Ferrari et al, Geophys. Res. 85 (1980) 3939. 

[17] J. O. Dickey, P. L. Bender, J. E. Faller, X. X. Newhall, 
R. L. Ricklefs, J. G. Ries, P. J. Shelus, C. Viellet, A. 
L. Whipple, J. R. Wiant, J. G. Williams, C. F. Yoder, 



Lunar laser ranging: A continuing legacy of the Apollo 
program, Science 265 (1994) 482. 

J. C. Ries , R. J. Eanes, C. K. Shum, M. M. Watkins, 
Geophys. Res. Lett. 19 (1992) 529. 

P. Dunn , M. Torrence, R. Kolenkiewicz, D. Smith, Earth 
scale defined by modern satellite ranging observations, 
Geophy. Res. Lett. 26 (1999) 1489. 

Y. Sofue and V. Rubin, Rotation curves of spiral galaxies, 

Ann. Rev. Astron. Astrophy. 39(2001)137. 

A. F. Zakharov et al., Solar system constraints on R n 

gravity, Phys. Rev. D74 (2006) 107101. 

A. L. Erickcek, T. L. Smith and M. Kamionkowski, ,Solar 

system tests do rule out 1/R gravity, Phys. Rev. D 74 

(2006) 121501 [arXiv:astro-ph/0610483l . 

T. Chiba, 1/R gravity and scalar-tensor gravity, Phys. 
Lett. B 575 (2003) 1 [arXiv:astro-ph /0307338] . 
K. Kainulainen, J. Piilonen, V. Reijonen and D. Sun- 
hede, Spherically symmetric spacetimes in f(R) gravity 
theories, Phys. Rev. D 76 (2007) 024020 [arXiv: 0704. 27291 

[gr-qc]]. 

G. J. Olmo, Limit to general relativity in f(R) the- 
ories of gravity, P hys. Rev. D 75, 023511 (2007) 
[arXiv:gr-qc/0612047j . 

S. DeDeo and D. Psaltis, Stable, Accelerating Universes 
in Modified Gravity, arXiv:0712.3939 [astro-ph]. 
A. Rajaraman, Newtonian gravity in theories with in- 
verse powers of R, arXiv:astro-ph/0311160 
T. Multamaki and I. Vilja, Spherically symmetric 
solutions of modified field equations in f(R) theo- 
ries of gravity, Phys Rev. D 74 (2006) 064022 
[arXiv:astro-ph/0606373] . 

V. Faraoni, Solar system experiments do not yet veto 
modifi ed gravity mod els, Phys. Rev. D 74 (2006) 023529 
[arXiv:gr-qc/0607016] . 

M. L. Ruggiero and L. Iorio, Solar System planetary or- 
bital motions and f(R) Theories of Gravity, JCAP 0701 

(2007) 010 [arXiv:gr-qc/0607093j . 

G. Allemandi, M. Francaviglia, M. L. Ruggiero and 
A. Tartaglia, Post-Newtonian parameters from alterna- 
tive theories of gravity, Gen. Rel. Grav. 37 (2005) 1891 
[arXiv:gr-qc/0506123] . 

Q. Exirifard, Phenomenological covariant approach to 
gravity, larXiv:0808T962l [gr-qc] . 

R. Saffari and S. Rahvar, f(R) Gravity: From the Pioneer- 
Anomaly to the Cosmic Acceleration, Phys. Rev. D 77 

(2008) 104028 |arXiv:0708. 14821 [astro-ph]]. 



